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1.  Introduction 


1.1  Objective  and  Scope  of  Effort 

This  work  was  performed  in  response  to  topic  AF93-015,  “Structmal  Integrity  of 
Intelligent  Materials  and  Structures.”  Intelligent  materials  open  new  avenues  to  im¬ 
prove  performance  reliability,  and  longevity  of  future  aerospace  vehicle  structures  by 
allowing  the  structural  materials  themselves  to  become  active  elements  for  multiple 
system  functions.  However,  the  application  of  these  materials  and  structures  has  been 
inhibited  because  the  effects  of  microstructural  interactions  between  intelligent  and 
host  material  elements  have  not  yet  been  well  characterized.  The  research  described 
here  was  to  develop  predictive  models  of  these  local  phenomena  to  enable  analytical 
assessment  of  their  effects  on  material  and  structural  integrity.  The  work  focused  on 
interfacial  loadings  and  propagation  of  interface  failures  in  intelligent  material  systems 
with  embedded  active  elements.  Detailed  analytical  and  finite-element  models  were 
developed  and  employed  to  assess  these  effects.  As  a  result,  some  key  characteristics 
governing  failure  mechanisms  in  these  systems  have  come  to  light. 

The  case  studies  encompassed  a  class  of  likely  application  scenarios  involving  em¬ 
bedment  of  piezoelectric  ceramic  wafers  in  conventional  polymer-matrix  composites. 
Two  particular  configurations  were  examined.  The  first  had  a  piezoceramic  replacing 
part  of  one  composite  layer  (“intralayer  embedment”)  while  the  second  had  the  ac¬ 
tive  wafer  inserted  between  layers,  thus  causing  a  bulge  in  the  material  (“interlayer 
embedment”).  Stress  and  strain  concentrations  were  examined  in  detail  for  regions 
near  interfaces  between  the  embedded  and  host  materials.  Finite-element  models 
were  used  to  assess  interfacial  crack  growth  behavior  both  between  the  piezoelectric 
element  and  the  host  material,  and  between  layers  of  the  host  material. 

The  results  of  these  studies  indicate  that  while  mode  I  (interlaminar  tension) 
responses  to  activation  of  the  piezoelectric  wafer  may  contribute  to  the  initiation  of 
interface  cracks,  crack  propagation  is,  at  least  for  the  cases  examined,  governed  by 
mode  II  (interlaminar  shear)  behavior  within  the  host  composite  and  at  the  interface 
between  intelligent  and  host  materials. 

1.2  Summary  of  Work  Performed 

The  work  done  tmder  this  project  has: 

1.  Developed  a  general,  two-dimensional  analytical  method  to  assess  stress  and 
strain  gradients  at  interfaces  within  intelligent  material  systems.  The  method 
was  implemented  in  a  computer  program  was  written  and  applied  successfully. 

2.  Developed  detailed  finite-element  models  of  composite  materials  incorporating 
embedded  intelligent  material  elements. 


3.  Implemented  the  virtual  crack  extension  method,  based  on  NASTRAN  finite 
element  models,  for  assessing  crack  growth  by  direct  computation  of  energy 
release  rates. 

4.  Examined  the  relative  severity  of  stress  and  strain  concentrations  subject  to  the 
influence  of  externally  applied  loads  as  well  as  induced-strain  activation  of  the 
intelligent  material  element 

5.  Assessed  the  influence  of  composite  lamination  angle  on  the  severity  of  local 
stress  gradients. 

6.  Assessed  factors  contributing  to  the  onset  and  propagation  of  interfacial  cracks 
in  intelligent  material  systems. 

7.  Considered  some  requirements  and  approaches  for  experimental  characterizar 
tion  of  interfacial  crack  propagation  in  intelligent  material  systems. 

Sections  2  through  5  discuss  the  four  tasks  that  were  defined  in  the  proposal, 
the  work  that  was  done,  and  the  results  that  were  obtained  for  some  representative 
configurations.  Section  6  presents  conclusions. 


2.  Task  1  -  Interface  Models 

Analytical  models  were  developed  to  allow  rapid  assessment  of  the  severity  of  local 
strain  gradients  caused  by  embedment  of  intelligent  material  elements  in  composite 
host  materials,  and  to  guide  development  of  finite-element  models  used  subsequently 
to  assess  fracture  behavior  of  these  systems.  The  analytical  models  were  based  on 
variable-order  power-series  expansions  of  structural  displacements  in  two  dimensions. 
High-order  models  were  employed  to  predict  strain  gradients  at  interfaces  between 
dissimilar  materials  in  the  intelligent  material  systems. 

A  variational  formulation  was  used  to  derive  the  governing  equations  for  the  an¬ 
alytical  model.  The  resulting  equations  for  the  unknown  coefficients  were  straight¬ 
forwardly  organized  into  a  linear  algebraic  system  and  solved  numerically.  Displace¬ 
ments,  strains  and  stresses  were  then  computed  from  the  power-series  expansions  and 
their  derivatives. 

The  variational  condition  is 


S{U  -  IV)  =  0  (1) 

where  U  and  W  are  strain  energy  and  work  terms,  respectively.  The  strain  energy 
is  expressed  in  terms  of  appropriate  derivatives  of  the  assumed  displacement  fields. 
The  work  term  involves  integrals  of  the  applied  surface  tractions  and  corresponding 
surface  displacements. 


U  =  (2) 

W  =  jj'^ddS 

Differentiating  the  above  integrals  with  respect  to  the  unknown  series  coefficients 
results  in  a  governing  equation  for  each  coefficient.  These  governing  equations  are 
assembled  into  an  algebraic  system  of  equations  and  solved  numerically. 

The  problem  addressed  under  this  project  considered  a  number  of  rectilinear  re¬ 
gions  adjoined  at  their  boundaries.  The  approach  was  to  divide  the  composite  mate¬ 
rial  system  into  such  regions  near  the  interfaces  of  interest.  Each  region  comprised  a 
single  material,  and  each  region  retained  its  own  local  displacement  expansion.  The 
regions  form  a  rectilinear  grid,  or  mesh,  wherein  each  individual  region  was  identi¬ 
fied  by  two  coordinate  indices  (I-J).  The  I-J  topology,  shown  in  Figure  1,  admitted 
modeling  of  both  embedded  and  siirface-attached  active-material  elements. 

Since  the  series  expansions  employed  are  continuous  beyond  first  derivatives,  they 
are  useful  only  within  a  single-material  region.  Thus,  matching  conditions  were  de¬ 
veloped  to  link  solutions  along  region  boimdaries  representing  interfaces  between  dis¬ 
similar  materials.  These  matching  conditions  were  implemented  using  Lagrange  mul- 


Local-Region  Reference  Frame 
and  Dimensions 


Figure  1:  I-J  Analytical-Model  Topology  (*  -  t,  j  indices  in  parentheses) 


tipliers  to  preserve  simplicity  and  numerical  precision  in  their  implementation.  The 
Lagrange  multiplier  method  also  3delds  exact  satisfaction  of  the  matching  conditions. 

The  two-dimensional  model  considered  surface  tractions  in  the  normal  and  tangen¬ 
tial  directions,  allowing  specification  of  all  admissible  boimdary-stress  components. 
Displacement  boimdary  conditions  were  limited  to  geometric  symmetry  conditions  on 
chosen  model  boundaries. 

The  following  sections  describe  details  of  the  analytical  model  formulation  and 
implementation. 


2.1  Analytical  Model  Formulation 

The  analytical  interface  model  developed  xmder  this  project  was  based  on  a  power- 
series  expansion  of  displacement  fields  in  two  dimensions.  Displacements  (u,w)  in 
the  x2-plane  were  assumed  to  be  of  the  form 

u{x,z)  =  aoo-l-Omo(®/0’"  +  aon(^/h)’*-»-amn(x/f)”‘(z/h)’* 

(3) 

w(x,z)  =  boo  +  bMer  +  boniz/hr-^b^nix/eri^^^^ 

The  quantities  £  and  h  axe  characteristic  spatial  lengths  in  the  x-  and  z-directions, 
respectively,  and  the  coefficients  a  and  b  are,  to  be  determined.  Equations  (2)  assume 
the  summation  convention  over  repeated  indices:  m  =  1 . . .  M,  n  —  1 . . .  N. 

Using  (3),  strains  linear  and  boundary  displacements  were  derived  for  a  chosen 
order  of  approximation  (A/,  IV).  The  linear,  total  strain  tensor  was  defined,  using 
contracted  notation,  as 


where  the  Cj  are  mechanical  strains,  and  the  ef  are  pre-strains  induced  by  thermal 
and/or  electrical  fields  (assumed  constant  within  each  region). 

Setting 

A  =  ^,  bmn  =  0  {m>  M  or  n>  N), 

displacements  and  their  spatial  derivatives  can  be  written 

u  =  aoo  +  amo(Aa:)"‘  +  aon(C-z)”  +  <imn(Ax)”‘(C2)” 
u„  =  A  [flio -I- ai„(Cr)"]  +  A(Tn -f- 1)  a(m+a)o(Ax)”* -I- a(,„+i)„(Ax)'"(Cr)”] 
u,z  =  C[<ioi +ami(Ax)"*] -t-C(n+ 1)  ao(n+i)(C2)”  +  am(n+i)(Ax)”‘(Cr)"] 
w  =  6oo  +  6mo(Ax)’”-H6on(C2)”  +  &mn(Aa:)”*(C2)’* 

W,x  —  A  [610 -f- 6ln(C^)”]  +  A(m -I- 1)  ^6(m+l)o(Ax)"‘ ■+•  6(m+l)n(Ax)”*(Cr)”] 

=  C(^)l  +  *ml(Ax)”*[-|-C(«-t-l)[6o(„+l)(C2)"  +  fem(n+l)(Ax)”*(C2)"] 

Taking  the  liberty  to  define  0°  =  1  allows  these  quantities  to  be  expressed  more 
compactly  as 

u  =  amn(Ax)'"(Cz)’* 
u,x  =  A(m -1- l)a(,n+i)n(Ax)”’(Cz)’* 
u,z  =  C(^+ l)am(n+i)(Ax)"‘(Cz)’* 
w  =  6mn(Ax)"*(Cz)’* 

VJ,x  =  A(m -H  l)6(m+l)n(Ax)”*(Cz)’* 

W,z  =  C(n+l)6m(n+l)(Ax)’"(C2)’* 


(m  =  0 . . .  M,  n  =  0 . . .  N).  (4) 


The  mechanical  strains  are  then 


75 


u„-c? 

Wyz-4 

U,z+W,x 


(5) 


where  contracted  notation  is  used  for  the  strain  components  and  the  quantities  ef  and 
ef  are  pre-strains  representing  a  combination  of  thermal  and/or  electrically  induced 
strains. 

The  analytical  model  assmned  a  transversely  isotropic  material  with  isotropy  in 
the  xz  plane.  Th\is,  the  strain  energy  within  a  region  was  expressed  as 

+  2Ci3(u,*  -4)  +  C'ssTs] 

Applied  tractions  were  also  expressed  as  power-series  expansions 

ai(£,z)  =  Sfc(C2)*'  ' 
fs(£,z)  = 

-r  M  .y  .J  }(j  =  0...M,  k  =  0...N}.  (7) 

az{x,h)  =  Tj{\xy 

n{x,h)  =  qj{\xy 

allowing  the  work  integral  to  be  expressed  as 

^  =  !Q[ci{£,z)'u{£,z)-{-f5{£,z)’w{i,z)]dz 

(8) 

+  So  [^3(2:,  h)  •  w(x,  h)  +  fs{x,  h)  •  u{x,  h)]  dx 

Inter-region  displacement  compatibility  and  symmetric  boundary  conditions  were 
enforced  by  means  of  Lagrange  multipliers.  The  compatibility  conditions  enforced 


=  u{£,zyi 
uix,  =  «(x, 

(9) 

u;(0, =  w{£,zy^ 

w{x,  =  w{x,  0)’^ 

The  boundary  conditions  enforced  were 

u(0,z)^J'  =  0 

w,.(0,.)‘'  =  0 

tx„(x,0)*L  =  0 
w{x,  0)*^  =  0 

The  resulting  Lagrangian  functions  xised  to  augment  the  energy  functional  were,  for 
the  matching  conditions, 


0 


-  u{e,  z)«]  = 

[tt;(x,  -  ty(x,  0)*^]  = 


and,  for  the  boundary  conditions, 

pi^u(0,2)^-'  =  0 

=  0 

(12) 

g'^u,z(x,oy^  =  0 
9wx^ixy0y^  =  U 

where  the  g^r,  are  unknown  Lagrange  multipliers. 

Expressions  (6),  (8),  (11)  and  (12)  were  assembled  into  the  energy  functional  (1) 
using  the  displacement  and  boundary-force  expansions  (3)  and  (7).  The  resulting 
integrals  were  evaluated  analytically,  and  differentiated  with  respect  to  each  of  the 
unknown  coefficients  amm  6mn  and  to  minimize  the  energy  functional. 

Analytical  integrals  of  the  energy  functional  constituents  were  easily  derived  from 
the  polynomial  forms  of  the  assumed  strain  fields  and  applied  boundary  forces.  Vol¬ 
ume  integrals  involve  terms  of  the  form 


r  f\xx)^+^{CzT'^^dxdz 

Jo  Jo 


for  products  of  coefficients  Omn  or  bmn  with  coefficients  ajk  or  bjk-  These  terms 
integrate  to 

r  f‘(xxr*Hczr*^dxdz  = 

Jo  Jo  ^  AC(m-l- j -H  l)(n-l-fc-l- 1) 

for  the  product  of  the  and  coefficients. 

The  work  terms  (surface  integrals)  involve  integrals  like 


f  «fc53amn(C^)"'^*  dz. 
L  m 


with  the  result 


I  S  amn(C‘2)’*'^^  dz  =  Y, 

•'®  L  m  m 


afcttmn 

C(n-l-fc+  1)' 


Differentiating  the  above  integrals  with  respect  to  the  unknown  coefficients  results 
in  the  following  nontrivial  equations  for  each  region. 


^  f _ ^ Cl 

[c(m  +  i  -  l)(n  + A:  + 1)  C(n+1) 


d  1  yc 
dbmn  2  Jo  Jo 


bik  — 


A(m  +  j  +  l)(n  +  &  —  1)  ^  A(7n  +  1) 


r  I  cs 

' (m  +  i)(n  +  A:)  ({n  +  1) 


a 

9&mn 


Jo  /  -Ci)K z-4dxdz  = 

r  [  Cl 

L(”^  + 

)  r  1  /C 

“  o  /  /  C55(«,z+t«,*)dar(f2  = 

Tin  1.2  2o  JO 


a  1  /■'*  /C 

3 -  n  C55{u,^+W,:,)dxdz  = 

vdfufi  ^  »/o  Jq 

^nk  jn 

^  A(m  +  i  +  l)(n  +  A;  -  1)^^*  (m  +  j)(n  +  A:) 


1  rh  ft 

2  Jo  Jo  ^55(w,z+w,*)dx(i^  = 


—  fi  r  f'c 

dbmn  [2  Jo  Jo 
^  mk 


:ajk  + 


(m  +  j){n  +  k)  ^  X{m  +  j-l){n  +  k  +  l) 


_ _ 

C(n  +  A:  +  1) 

tk 

C(n  +  fc  +  l) 

Tk 

X{m  +  fc  + 1) 


9 


L|f^||iM«r|n||^n^ 

2.1  Analytical  Model  Formulation 


d 

dbmn 


dx 


Qk 

A(J71  +  /k  +  1) 


The  compatibility  constraint  relations  (11)  take  the  following  form: 


sS"‘’(aS'‘’-E«<in)  =  0 

=  0 

sSt'Ki>M'’  -'ZntfL)  =  0 

The  boundary  constraint  relations  (12)  take  the  form: 


(14) 


9utrfimX 

/iO  }iX 

These  expressions  are  straightforwardly  differentiated  with  respect  to  Omm  i>mn  and 
g^ri  coeflScients. 

The  above  differentiated  expressions  comprise  a  set  of  2{M  +  1)(JV'  + 1)  equations 
for  the  displacement-expansion  coefficients  Omn  and  6mn  in  each  region.  In  addition, 
each  region  is  associated  with  2(M-l-iV-t-2)  Lagrange  multipliers  to  enforce  boimdary 
and  matching  conditions.  These  equations  are  assembled  into  a  symmetric,  linear 
algebraic  system  of  equations  for  the  entire  structure  to  be  modeled,  and  the  system 
of  equations  is  solved  subject  to  prescribed  boundary  forces  and/or  imposed  internal 
strains  (ef). 

The  order  of  the  displacement-expansion  approximation  is  set  by  choosing  values 
M  and  N  for  each  region,  with  the  restrictions  that  Mij  be  the  same  for  all  regions 
with  the  same  i-index,  and  that  Nij  be  the  same  for  all  regions  with  the  same  j-index. 
In  practice,  M  and  N  values  were  set  the  same  for  all  modeled  regions,  although  M 
was  not,  in  general,  set  to  equal  N. 

One  sho\ild  note  that  the  parameters  M  and  N  control  the  fidelity  of  the  solu¬ 
tion  provided  by  the  analytical  model.  In  general,  these  solutions  will  be  subject 
to  tnmcation  errors  unless  very  high-order  approximations  are  used.  In  the  present 
effort,  convergence  behavior  of  the  anal}rtical  model  was  studied,  and  resulted  in  the 
selection  of  Af  =  9,  =  4  for  the  material  configurations  considered  (see  §3.1). 


0 

0 

0 

0 
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2.2  Analytical  Model  Implementation 

The  analytical  model  described  in  the  previous  section  was  implemented  in  a  C- 
language  computer  program.  The  program  is  written  in  standard,  ANSI-C,  and  may 
be  nm  on  any  computer  having  an  ANSI-C  compiler.  For  this  project,  the  analytical 
model  was  implemented  on  an  Apple  Macintosh  Ilci  computer. 

Two  principal  advantages  were  found  in  using  the  analytical  model.  First,  since 
the  model  addresses  only  a  specific  problem,  it  executes  much  faster  than  models  of 
equivalent  fidelity  implemented  in  commercial  finite-element  programs.  This  makes 
the  analytical  model  practical  for  assessing  the  effects  of  configuration  variables  such 
as  material-layer  thicknesses  and  layup  angles.  Second,  the  model  converges  rapidly 
on  smooth  descriptions  of  the  high  strain  gradients  occurring  at  material  interfaces. 
This  allows  the  model  to  be  used,  as  it  was  in  the  present  study,  to  guide  the  dis¬ 
cretization  of  more  general,  finite-element  models  for  detailed  stress  and  fracture 
predictions  (see  §3.1.4). 

The  analytical  model  implemented  here  could  be  improved  for  faster,  more  gen¬ 
eral  operation.  For  example,  the  Lagrange  multiplier  implementation  of  boundary 
conditions  could  be  simplified  by  eliminating  prescribed  boundary  unknowns  directly 
before  assembly  of  the  governing  equations.  Also,  more  general  boimdary  and  inter¬ 
face  conditions  (e.g.,  internal  free  surfaces)  could  be  implemented  to  allow  modeling 
of  fracture  and  delamination  behavior.  Finally,  the  power-series  representation  of  dis¬ 
placement  fields  allows  analytical  integration  of  nonlinear  strains  to  be  implemented 
with  relative  ease.  This  would  allow  the  model  to  treat  geometrically  nonlinear  be¬ 
havior  as  well. 


Table  1:  Analytical  Model  Material  Properties 


Property 

Value 

SI  Units 

English  Units 

Composite  Moduli 

E\\ 

127  GPa 

18.4  Msi 

E33 

10  GPa 

1.5  Msi 

G31 

5.8  GPa 

0.85  Msi 

Ui3 

0.29 

0.29 

Piezoelectric 

Moduli 

fJii 

88  GPa 

12.8  Msi 

E33 

66  GPa 

9.6  Msi 

G31 

29  GPa 

4.2  Msi 

Vl3 

0.31 

0.31 

Induced  Strain: 

di3 

-122  pC/N 

dss 

285pC/N 

3.  Task  2  —  Models  of  Embedded  Elements 

The  analytical  model  described  in  the  previous  section  was  employed  along  with 
detailed  finite-element  models  to  assess  structmal  responses  and  fracture  behavior 
near  interfaces  between  intelligent  and  host  composite  material  elements.  The  models 
employed  in  this  investigation,  and  results  obtained  from  them,  are  described  in  this 
section.  The  results  presented  consist  of  predictions  of  interfacial  strain  concentrations 
for  various  host-material  orientations,  assessment  of  discretization  requirements  for 
finite-element  models  of  the  interface  regions,  and  finite-element  predictions  of  local 
strain  gradients  near  material  interfaces. 

3.1  Anal3rtical  Models 

3.1.1  Description 

The  analytical  model  was  applied  to  the  configuration  shown  in  Figure  2  to  assess  local 
strain  distributions.  This  model  represents  a  piezoelectric  wafer  embedded  within  a 
layer  of  a  thin,  high-strength,  graphite-epoxy  laminate.  The  properties  employed  for 
each  material  are  listed  in  Table  1. 

The  model  was  loaded  by  a  combination  of  piezoelectric  induced  strains  and  an 


Figure  2:  Analytical  Interface-Stress  Model  Configuration 

axial  stress  of  100  MPa  (14.5  ksi).  The  piezoelectric  strains  corresponded  to  an  applied 
electric  field  of  400  kV/m  of  across  the  0.254  mm-thick  wafer,  resulting  in  applied 
strains  e?  =  —48  •  10“®  and  C3  =  112  •  10“®. 

The  analytical  models  were  exercised  for  two  purposes.  First,  the  analytical  mod¬ 
els  were  \ised  to  determine  the  effect  of  lamination  angle  on  the  severity  of  resulting 
strain  gradients.  This  analysis  was  xmdertaken  to  study  design  modifications  that 
might  result  in  less-severe  strain  concentrations  at  material  interfaces.  Second,  strain 
gradients  were  predicted  to  determine  smtable  discretizations  for  subsequent  finite- 
element  analyses.  This  was  imdertaken  ensure  that  finite-element  models  employed  to 
assess  crack-growth  behavior  in  the  intelligent  material  system  would  provide  accurate 
results. 

3.1.2  Interface  Strain  Field  Predictions 

The  configuration  shown  in  Figvure  2  was  analyzed  for  two  cases:  an  applied  axial 
stress  loading  of  100  MPa,  and  applied  piezoelectric  strains  €?  =  -48  •  10”®  and 
ef  =  112  •  10”®.  The  results  of  these  analyse  are  shown  in  Figures  3  to  8.  Figures  3 
to  5  show  strain  distributions  along  the  length  of  the  interlaminar  interface  at  the 
top  of  the  piezoelectric  wafer  (z  =  0.127  mm)  as  a  function  of  distance  from  the 
symmetry  plane  imposed  at  x  =  0.  Figures  6  to  8  show  strain  distributions  through 
the  thickness  of  the  model  at  the  interface  between  the  embedded  piezoelectric  wafer 
and  the  host  material  (x  —  0.025  m)  and  from  the  midplane  (z  =  0)  to  the  upper 
(free)  surface.  In  these  figures,  the  legend  “Force  Loading”  refers  to  application  of 
the  100  MPa  axial  stress,  and  the  legend  “Strain  Loading”  refers  to  the  case  where 
piezoelectric  strains  are  prescribed.  The  ‘‘upper  face”  designation  refers  to  strains 
in  the  composite  host  material  just  above  the  interface  (z  =  0.127+  mm),  while 
the  “lower  face”  designation  refers  to  strain  distributions  just  below  the  interlaminar 
interface  (z  =  0.127—  nun).  The  “+X”  and  “-X”  designations  refer  to  locations  either 
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Table  2:  Interfadal  Strain  Concentrations 


Strain  Measure  (%) 

100  MPa  Loading 

eP  Loading 

0.012 

0.001 

0.029 

0.010 

ItT^I 

0.066 

0.0092 

side  of  the  interface  at  (x  =  0.025  m),  with  the  “-X”  referring  to  the  side  containing 
the  piezoelectric  wafer. 

Examination  of  Figures  3  to  8  reveals  some  important  aspects  of  the  material 
system’s  elastic  behavior.  First  and  foremost,  strain  concentrations  at  the  interfaces 
are  influenced  far  more  by  the  application  of  the  100  MPa  mechanical  stress  than  by 
the  applied  piezoelectric  strains.  Naturally,  the  magnitudes  of  these  concentrations 
will  depend  directly  on  the  magnitude  of  the  applied  loading.  However,  the  100  MPa 
applied  stress  represents  roughly  10%  of  the  composite  material’s  ultimate  strength, 
while  the  piezoelectric  applied  strains  represent  roughly  50%  of  that  material’s  actu¬ 
ation  capability.  Table  2  lists  peak  strain  concentrations  for  the  two  loading  cases. 
Examination  of  this  table  shows  that  the  integrity  of  the  intelligent  material  sys¬ 
tem  will  be  most  affected  by  the  external  loading  environment.  Only  in  cases  where 
external  loads  are  small  will  induced-strain  effects  cause  the  worst  interfacial  stresses. 

Table  2  shows  that  shear  strain  concentrations  are  the  most  severe  under  applied 
in-plane  loadings,  but  through-thickness  strain  concentrations  are  slightly  more  severe 
when  only  induced-strain  loading  is  present.  Prom  these  results,  one  might  expect 
the  action  of  the  intelligent  material  to  have  an  appreciable  effect  on  mode  I  (tension) 
interlaminar  crack  initiation,  but  subsequent  crack  growth  would  be  heavily  influenced 
by  mode  II  (shear)  behavior. 

Another  observation  to  be  drawn  from  the  model  predictions  is  that  the  char¬ 
acteristic  lengths  of  the  strain  concentrations  are  very  small,  typically  of  the  order 
of  1-1.5  mm.  (4-6  layer  thicknesses).  Also,  both  through-thickness  (€3)  and  shear 
(75)  strain  gradients  exhibit  reversals  with  increased  distance  from  the  interface.  The 
magnitude  of  the  through-thickness  strain  reversal  is  roughly  0.002%  (one-tenth  of 
the  strain  concentration  magnitude),  and  the  magnitude  of  the  shear  strain  reversal 
is  roughly  0.005%  (less  than  one-tenth  of  the  strain  concentration  magnitude). 


Through-Thickhess  Strain  (mechanical) 


Distance  from  Centerline  (m) 

Figure  3:  Axial  Strains  along  Interlaminar  Interface 


Distance  from  Centerline  (m) 

Figure  4:  Through-Thickness  Strains  along  Interlaminar  Interface 


Distance  from  Laminate  Midplane  (m) 


3.1  AnsJytical  Models 


Distance  from  Centerline  (m) 

Figure  5:  Shear  Strains  along  Interlaminar  Interface 


Figure  6:  Interface  Axial  Strains  through  Laminate  Thickness 


Distance  froin  Laminate  Midplane  (m)  1  Distance  from  Laminate  Midplane  (m) 


3.  TASK  2 -MODELS  OF  EMBEDDED  ELEMENTS 


Through-Thickness  Strain  (mechanical) 

e  7:  Interface  Through-Thickness  Strains  through  Laminate  Thickness 


Figure  8:  Interface  Shear  Strains  through  Laminate  Thickness 
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Composite  Layup  Angle  (deg.) 


Figure  9;  Strain  Concentration  vs.  Lamination  Angle  -  Applied  Stress  Loading 


3.1.3  Influence  of  Host-Material  Orientation 

Subsequent  analyses  employing  the  analytical  model  were  run  to  assess  the  effects  of 
composite  lamination  angle  on  the  severity  of  stress  gradients  in  the  intelligent  mate¬ 
rial  system.  The  purpose  of  these  analyses  was  to  recommend  design  configurations 
that  might  minimize  interfacial  stress  concentrations. 

Figures  9  and  10  show  maximum  through-thickness  (normal)  and  shear  inter¬ 
facial  stress  concentrations  versus  lamination  angle  of  the  host  composite  material 
(note  that  the  composite  laminate  is  balanced  with  [±d\  plies  to  accommodate  the 
two-dimensional  assumptions  of  the  analytical  model).  These  figures  show  that  both 
normal  and  shear  stress  concentrations  grow  under  applied  axial  loadings  as  the  lam¬ 
inate  axial  stiffiiess  drops.  However,  normal  stress  concentrations  are  relatively  in¬ 
sensitive  to  lamination  angle  under  activation  of  the  intelligent  material  wafer.  Shear 
stress  concentrations  actually  reduce  with  increasing  off-axis  laminate  angle  under 
induced-strain  loading. 

The  effect  of  the  variations  shown  in  Figures  9  and  10  is  to  recommend  axially 
dominated  laminates  for  axial  loadings.  As  noted  in  §3.1.2,  interface  normal  stresses 
are  expected  to  dominate  for  initiation  of  interlaminar  cracks  under  applied,  in-plane 
loadings.  Aligning  the  composite-material  layers  with  this  loading  provides  minimtim 
normal-  and  shear-stress  concentrations;  Furthermore,  no  appreciable  reduction  in 
normal-stress  concentration  can  be  gained  by  off-axis  lamination  even  in  the  case 
where  only  induced-strain  loadings  are  present,  as  seen  in  Figure  10.  However,  if  the 
material  loading  is  to  be  dominated  by  strain  activation  of  the  intelligent  material,  re¬ 
ductions  in  shear-stress  concentrations  could  be  achieved  by  increasing  the  composite 


Composite  Layup  Angle  (deg.) 

Figure  10:  Strain  dloncentration  vs.  Lamination  Angle  -  Induced  Strain  Loading 
lamination  angle. 

Another  instructive  way  to  view  these  results  is  with  respect  to  the  relative  stifhess 
of  the  composite  host  material  in  axial  and  through-thickness  shear  directioiis.  As 
the  lamination  angle  is  increased,  the  transverse  shear  stiSuess  C55  increases  relative 
to  the  axial  stiffness  Cn.  Figures  11  and  12  show  the  effect  of  this  stiffness  ratio. 
These  figures  indicate  that  use  of  stiffer  composite  matrix  materials  raise  normal-  and 
shear-stress  concentrations  xmder  applied  external  loadings,  and  lower  the  shear-stress 
concentration  under  pure  induced-strain  loading.  Accordingly,  if  significant  external 
loadings  are  to  be  applied  to  the  intelligent  material  system,  interfacial  stresses  are 
reduced  as  compliance  of  the  host-composite  matrix  is  increased. 

3.1.4  Requirements  for  Model  Discretization 

The  analyses  presented  in  §3.1.2  were  also  employed  to  determine  suitable  discretiza¬ 
tions  for  finite-element  models  used  in  crack-growth  studies.  Two  discretization  re¬ 
quirements  were  enforced.  First,  the  finite-element  model  was  required  to  span  a 
length  sufBcient  to  capture  the  entirety  of  the  strain  gradients  either  side  of  the  ma¬ 
terial  interface.  Examination  of  Figures  3  to  5  shows  this  length  to  be  ±  3  to  5  mm 
from  the  interface. 

Second,  the  finite-element  model  was  requued  to  have  sufficiently  small  elements 
to  enable  capture  of  the  strain-gradient  magnitudes  at  the  material  inter&ce.  For 
linear-displacement  (constant  strain)  elements,  a  minimiiTn  of  two  elements  are  re¬ 
quired  to  capture  a  linear  strain  gradient.  lb  match  the  8^-order  gradient  predicted 
by  the  analytical  model,  a  minimum  of  nine,  constant-strain  elements  would  be  re- 
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Figure  13:  Intralayer  embedded  FEM  model 

quired.  Capturing  the  strain  gradient  predicted  by  the  analytical  model  thus  requires 
a  minimum  of  nine  elements  over  the  first  millimeter  on  either  side  of  the  interface. 
Maximiim  element  dimension  according  to  this  criterion  is  roughly  0.1  mm,  or  about 
40%  of  a  layer  thickness. 

Through-thickness  gradients  must  also  be  captured  in  the  finite-element  model. 
To  capture  the  third-order  strain  field  of  the  analytical  model  requires  a  minimum  of 
four  constant-strain  elements  through  the  thickness  of  each  layer.  This  translates  into 
a  maximum  element  dimension  of  roughly  0.032  mm  (0.00125  in),  or  about  one-third 
of  the  minimtun  required  axial  dimension. 

An  element  dimension  of  0.25  •  10~^  inch  (0.006  mm)  was  chosen  for  the  finite- 
element  models  used  in  this  study.  This  dimension  exceeded  the  discretization  re¬ 
quirement  established  by  the  analytical  model  by  over  a  factor  of  two,  as  the  layer 
thickness  employed  in  the  finite-element  model  was  half  that  modeled  analytically. 

3.2  Finite  Element  Models 
3.2.1  Intralayer  Embedment 
Description 

The  finite  element  model  used  the  configuration  shown  in  Figure  13.  This  model,  like 
the  analytical  model,  represents  a  piezoelectric  wafer  embedded  within  a  graphite- 
epoxy  laminate.  For  comparison  purposes,  the  material  properties  used  were  the 
same  orthotropic  properties  as  those  \ised  for  the  analytical  model.  All  layers  of  the 
griqjhite-epooQr  composite  had  a  fiber  lamination  angle  of  zero.  Tb  reduce  computer 
tion,  the  length  of  the  specimen  was  smaller  than  the  analjrtical  model. 


The  entire  model  was  constructed  of  plate  elements  6.35  fim  (0.25  •  10”®  in.)  on 
a  side  and  1  mm  (0.0394  in.)  thick.  Model  symmetry  was  used  so  only  the  upper 
right  quadrant  of  the  system  was  actually  modeled.  Coincident  nodes  were  used  along 
the  vertical  interface  between  the  piezoelectric  and  the  adjacent  horizontal  layer  of 
composite.  These  nodes  were  left  unattached  for  all  runs  of  the  model  because  no 
bonding  would  be  expected  along  this  interface.  Coincident  nodes  were  also  used 
along  the  horizontal  interface  between  the  top  of  the  piezoelectric  and  the  adjacent 
composite,  and  between  adjacent  layers  of  the  composite  at  the  same  height  as  the 
piezoelectric.  All  horizontal  coincident  nodes  were  attached  for  the  initial  runs.  The 
nodes  between  adjacent  composite  layers  were  then  released  in  subsequent  runs  to 
examine  crack  growth. 

An  axial  stress  of  1.75  MPa  (254  psi)  was  applied  to  right-hand  edge  of  the  model. 
Thermal  loadings  were  used  to  simulate  induced  piezoelectric  strains.  A  AT  value 
of  100  was  used  for  all  cases,  along  with  horizontal  and  vertical  thermal  expansion 
coefficients  of  -9.606E-7  and  2.2441E-6  (K)”^,  respectively.  This  allowed  a  direct 
correspondence  between  temperature  and  applied  voltage.  The  opposing  signs  of  the 
coefficients  provided  for  piezoelectric  expansion  along  one  axis  while  it  contracted 
along  the  other.  The  applied  piezoelectric  strain  for  all  cases  was  c?  =  —9.6  •  10”® 
and  cf  =  22.4  •  10”®. 

The  final  model  size  was  determined  by  building  a  long  model  and  examining 
the  locations  of  essentially  zero  shear  strain  given  the  combination  of  the  applied 
loads  above.  Figure  14  is  a  gray-scale  image  of  the  final  model  showing  that  the 
shear  strains  eflFectively  go  to  zero  within  three  thicknesses  of  the  piezoelectric  wafer. 
Also,  the  analytical  model  predicted  that  a  length  to  thickness  ratio  of  7.9  for  the 
piezoelectric  wafer  would  be  necessary  for  complete  stress  analysis.  The  final  FE 
model  has  a  length  to  thickness  ratio  of  10. 


Results  -  Stress  Gradients/Deformations 

Once  final  model  size  was  determined,  eleven  runs  were  performed  in  which  coincident 
nodes  along  the  composite  layer  interface  were  released  one  by  one.  An  initial  nm 
was  made  with  no  horizontal  nodes  released.  The  results  of  this  nm  supplied  dis¬ 
placements  and  node  restraint  forces  for  the  nm  in  which  the  first  coincident  nodes 
were  released.  Figures  15  and  16  show  the  nodal  displacements  given  by  these  first 
two  runs.  Work  terms  were  calculated  from  the  new  displacements  of  the  second  run. 
This  was  repeated  nine  more  times.  Figure  17  shows  the  displacement  field  after  four 
nodes  were  released.  As  expected,  the  percent  change  from  one  release  to  the  next 
for  the  mode  I  (normal-tension)  displacements  have  drastically  diminished  while  the 
percent  change  of  the  mode  II  (shear)  displacements  have  not  been  much  altered. 
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Figure  15:  Nodal  displacements  with  no  horizontal  nodes  released 
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Figure  17:  Nodal  displacements  with  four  horizontal  nodes  released 


Figure  18:  Interlayer  embedded  FEM  model 

3.2.2  Interlayer  Embedment 
Description 

The  model  used  in  the  previous  intralayer  embedment  was  modified  to  represent  a 
piezoelectric  wafer  embedded  between  two  layers  of  graphite-epoxy  laminate.  Fig¬ 
ure  18  shows  the  modified  model.  Coincident  nodes  were  placed  along  the  horizontal 
interface  along  the  top  of  the  piezoelectric  and  the  adjacent  composite.  These  nodes 
were  attached  for  initial  runs  and  then  released  in  subsequent  runs  to  examine  crack 
growth  between  the  piezoelectric  and  composite.  The  nodes  along  the  horizontal 
plane  of  symmetry  between  composite  layers  were  also  held  and  then  released  to 
study  crack  growth  between  composite  layers. 

The  previous  intralayer  configuration  used  an  axial  loading  of  1.75  MPa.  This 
was  done  to  try  to  equilibrate  the  effects  of  both  the  axial  load  and  the  induced 
piezoelectric  strain.  After  examination  of  the  results,  it  was  determined  that  the  axial 
load  was  the  prime  contributor  to  crack  growth  after  the  first  node  was  released.  It 
was  therefore  decided  to  increase  the  magnitude  of  the  axial  load  for  the  interlayer 
configuration  100-  fold  to  175  MPa  (25.4  ksi).  This  was  applied  in  separate  runs  for 
both  tension  and  compression.  The  induced  piezoelectric  strains  were  not  changed 
in  magnitude  from  those  used  in  the  intralayer  configuration,  although  the  sign  of 
AT  was  reversed  for  all  runs  under  axial  compression.  This  configuration  -  axial 
compression  loading  combined  with  negative  piezoelectric  strain  in  the  Z  direction  - 
gave  the  highest  nodal  restraint  forces. 


Results  -  Stress  Gradients/Deformations 
A  total  of  fifteen  runs  were  performed: 


•  five  nins  with  tensile  axial  load  in  which  single  point  constraints  for  the  nodes 
along  the  line  of  horizontal  symmetry  were  released  one  by  one 

•  five  runs  with  compressive  axial  forces  in  which  the  coincident  nodes  along  the 
piezoelectric-composite  layer  interface  were  released  one  by  one 

As  with  the  intralayer  configuration,  an  initial  run  was  made  for  both  the  tensile 
and  axial  compression  load  cases  with  no  nodes  released.  Then,  each  successive  run, 
tensile  as  well  as  compressive  load  cases,  built  on  the  rxm  before  it  to  detennine  work 
terms  (release  force  •  relative  displacement)  in  both  the  X  and  Z  directions  for  each 
node  released.  Figures  20  and  21  show  the  model  in  tensile  and  compressive  loadings, 
respectively. 

The  interesting  note  from  the  tensile  case  is  that  the  node  release  force  for  the  first 
node  between  the  composite  layers  was  positive  in  the  Z  direction.  This  indicates  that 
at  this  point,  where  crack  initiation  might  be  expected,  there  is  actually  a  compressive 
force  that  inhibits  mode  I  cracking.  However,  subsequent  nodes  to  the  right  of  this 
experience  a  tension  force  tending  to  cause  cracking. 

Note  also  that  the  piezoelectric  in  Figure  20  has  been  displaced  in  the  positive  X 
direction  even  though  the  induced  strain  for  this  case  was  in  the  positive  Z,  negative 
X  direction.  This  positive  displacement  can  be  attributed  entirely  to  the  large  applied 
axial  tensile  stress  that  overwhelms  the  induced  piezoelectric  strain.  Similarly,  the 
piezoelectric  in  Figure  21  has  been  displaced  in  the  negative  X  direction  even  though 
the  induced  strain  for  this  case  was  in  the  positive  X  direction. 

The  axial  compression  load  condition  was  the  more  interesting  as  the  node  release 
forces  for  mode  II  crack  growth  were  an  order  of  magnitude  higher  than  the  next 
closest  release  rates  of  any  configuration.  Figure  22  shows  the  displacement  field 
under  compressive  axial  load  after  four  nodes  were  released  between  the  piezoelectric 
and  composite  interface. 


Figure  19:  Nodal  displacements  with  no  nodes  released;  axial  load  in  tension 
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Figure  20:  Nodal  displacements  with  no  nodes  released;  axial  load  in  compression 


Figure  21:  Nodal  displacements  with  four  nodes  released;  axial  load  in  compression 


Figure  22;  Quantities  for  Calculating  Strain  Energy  Release  Rates 


4.  Task  3  —  Interface  Failure  Models 

The  finite-element  models  described  in  §3.2  were  used  to  assess  crack  propagation 
behavior  in  the  intelligent  material  system.  The  method  employed  for  this  assess¬ 
ment  computed  strain-energy  release  rates  for  growing  cracks  by  calculating  the  work 
done  by  crack-closure  forces  associated  with  nodal  displacement  constraints  along  a 
predefined  crack  path.  Energy  release  rates  were  investigated  for  two  cases;  one  in 
which  the  intelligent  material  element  (piezoelectric  wafer)  was  inserted  into  a  layer 
of  the  host  composite  material  (“intralayer”  embedment),  and  one  in  which  the  in¬ 
telligent  material  element  was  inserted  between  layers  of  host  material  (“interlayer” 
embedment).  In  the  former  case,  energy  release  rates  were  calculated  for  a  combined 
piezoelectric  actuation  and  applied  axial  tension  loads.  In  the  latter  case,  both  ten¬ 
sile  and  compressive  axial  loads  were  applied  in  conjimction  with  the  piezoelectric 
actuation  to  assess  crack  propagation  at  two  sites  in  the  material  S3rstem. 

4.1  Calculation  of  Strain-Energy  Release  Rates 

Strain-energy  release  rates  were  computed  from  the  finite-element  model  data  \ising  a 
node-release  method.  This  method  computed  the  work  released  by  releasing  a  nodal 
displacement  constraint  along  the  crack  path,  as  illustrated  in  Figure  22.  The  method 
involves  comparing  results  of  two  analysis  steps.  First,  constraint  forces  are  recorded 
when  the  crack  is  held  closed  by  a  constraint  applied  to  enforce  displacements  at 
coincident  nodes  along  the  crack  path  to  be  equal.  Second,  the  constraint  is  released 
and  the  resulting  relative  nodal  displacements  are  recorded.  The  work  released  by 
growing  the  crack  to  the  next  coincident-node  pair  along  the  cradc  path  is  then 
computed  as 


(16) 


I.  Iwc. 


The  node-release  method  allows  the  relative  contribution  of  mode  I  (normal  ten¬ 
sion)  and  mode  II  (shear)  stress  intensities  to  be  assessed.  The  first  work  term  on 
the  right-hand-side  of  eq.  (16)  represents  the  shear  (mode  II)  contribution,  while  the 
second  represents  the  tensile  (mode  I)  contribution.  The  strain-energy  release  rates 
are  thus 

G,  =  (17) 

Gn  =  i/xAu. 

For  very  long  cracks,  an  expression  derived  by  Wisnom^  can  be  employed  to  assess 
the  asymptotic,  mode  II  strain-energy  release  rate: 

r 

"  4Fu(h  - 1) 

where  (j\  is  the  applied  axial  stress,  h  is  the  total  thickness  of  the  laminate,  t  is 
the  thickness  of  the  material  region  between  the  growing  cracks  on  either  side  of 
the  midplane  of  the  laminate,  and  En  is  the  composite  axial  modulus.  For  the 
present  models,  which  represent  one  half  of  Wisnom’s  model  through  the  thickness, 
the  asymptotic  strain-energy  release  rate  is  one-half  that  given  in  eq.  (18),  and  takes 
a  value  of 

2 

Giis  =  1.875  •  10"®^  in  (19) 

•C'll 

=  4.763- 10”®^  m 
E\\ 

for  the  properties  and  dimensions  employed  here.  It  is  useful  to  examine  the  value 
of  Gii  relative  to  its  asymptotic  value  to  determine  how  far  a  delaminating  crack 
must  grow  before  the  delamination  propagates  catastrophically. 

4.2  Virtual  Crack  Extension 

4.2.1  Theory 

The  virtual  crack  extension  method^  may  be  used  to  compute  strain  energy  release 
rates  along  highly  stressed  interfaces  in  material  systems  such  as  those  studied  here. 

^Wisnom,  M.  R.,  “Delamination  in  Tapered  Unidirectional  Glass  Fibre-Epoxy  under  Static  Ten¬ 
sion  Loading,”  AIAA-91-1142-CP,  Prvc.  SSnd  Structures,  Structural  Dynamics  and  Materials  Con¬ 
ference,  April  1991,  p.  1162 

^HeUen,  T.  K.,  “On  the  Method  of  Virtual  Crack  Extensions,”  Int.  J.  Num.  Method.  Eng.,  Vol. 
9  (1975)  p.  187. 
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Virtual  Crack  Extension  Aa 


a  =  j^lt/(d  +  Aa)-(/(d-A»)l 
-j|^iV'(d  +  Aa)-V(d-A«)l 


Figure  23:  Virtual  Crack  Extension  Computation  of  Strain  Energy  Release  Rate 

Evaluation  of  these  release  rates  can  lead  to  predictions  of  the  likelihood  of  compro¬ 
mising  material  and  structural  integrity  by  propagating  interfacial  cracks  and  delam¬ 
inations  into  the  host  material. 

The  method  relies  on  a  relatively  good  description  of  the  stress  state  in  a  region 
somewhat  removed  from  a  hypothesized  crack  tip  to  compute  an  approximation  to 
the  instantaneous  strain  energy  release  rate  G.  If  this  rate  exceeds  the  material 
limit,  the  crack  is  grown  to  the  next  mesh  point  and  the  analysis  is  continued.  Some 
recent  research  has  shown  very  good  agreement  between  analytical  VCE  models  and 
experiments  for  delamination  in  certain  composite  constructions.® 

Practically  speaking,  virtual  crack  extension  must  be  implemented  in  the  context 
of  a  finite  element  model  with  good  detail  around  the  crack  tip,  as  illustrated  in 
Figure  23.  A  stress  analysis  is  performed,  resulting  in  a  set  of  nodal  displacements 
d.  A  ring  of  elements  is  selected  that  surrounds  the  crack  tip  but  is  several  elements 
away  from  it.  An  energy  release  rate  G  is  computed  by  applying  a  small  rigid-body 
displacement  to  the  nodes  on  the  inside  of  the  ring,  shown  as  a  horizontal  motion 
in  the  figure  for  mode  I  crack  growth.  Thus  if  r  is  a  vector  of  unit  motions  in  the 
specified  direction  along  the  inner  ring  of  nodes,  we  have 

G  =  ^  [U(d  +  Aar)-  Cf  (d  -  Aa  r)]  -  [V{d  -f  Aa  r)  -  V(d  -  Aa  r)]  (20) 

where  U  represents  the  strain  energy  developed  as  a  result  of  this  fictitious  movement 
^Wisnom,  op.  dt. 


i.2..  .Yktuai  Crack  Ex^osiqq 


and  V  represents  work  done  by  external  loads. 
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4.2.2  Development 

Assuming  no  loads  applied  to  the  ring  of  elements  and  hence  no  work  term  V,  the 
required  strain  energy  values  can  be  obtained  using  the  stifihess  matrix  for  the  ring 
of  elements: 

1 


Expanding, 


so  that 


£;(d  +  Aar)  =  ^(d  + Aorf  K(d  + Aar) 

A 

C/(d-Aar)  =  ^(d  -  Aar)’'K(d  -  Aar) 

A 


J7(d  +  Aa)  =  ^d^Kd  + Aad^Kr  +  ^(Aa)^r^Kr 
£/(d  -  Aa )  =  ^d^Kd  —  Aa  d^Kr  +  ^(Aa  )^r^Kr 

A  ^ 


2Aa 


(21) 


(22) 

(23) 


a  result  which  is  seen  to  be  independent  of  the  displacement  increment  Aa. 

Note  that  the  vector  of  nodal  point  forces  generated  by  the  elements  in  the  ring 
is  f  =  Kd,  so  that  the  expression  reduces  to 


G  =  rr 


(24) 


Recognizing  that  r  is  a  vector  of  unit  values  corresponding  to  translations  of  each  of 
the  two  rings  of  nodes  in  a  specified  direction  (horizontal  in  Figure  23),  we  see  that 
G  is  simply  the  sum  of  the  forces  applied  by  the  elements  in  the  ring  to  the  nodes  on 
the  inner  ring.  The  forces  fi:om  one  side  of  the  ring  tend  to  cancel  those  on  the  other 
side,  and  this  suggests  that  considerable  precision  may  be  required  in  performing  this 
calculation. 

The  foregoing  derivation  assinned  that  no  loads  were  applied  to  any  elements  in 
the  ring.  However,  as  explained  elsewhere,  piezoelectric  actuation  was  simulated  in 
the  finite  element  models  used  in  this  project  by  thermally  induced  strsdns.  Thus  there 
were  in  fact  loads  applied  to  those  elements  in  the  ring  that  represent  piezoelectric 
material.  In  order  to  account  for  this  one  must  replace  f  by  f + g  where  g  is  a  vector 
of  nodal  loads  due  to  the  fictitious  temperature  rise.  These  load  vectors  are  generated 
by  finite  element  codes  by  integration  of  the  thermal  stresses  of  the  form  Ea/\T  over 
the  area  of  the  element.  The  energy  release  rate  is  then 

G  =  (f  +  gfi 


(25) 


Note  that  this  term  must  be  divided  by  the  thickness  of  the  elements  used  to  represent 
the  “slice”  of  composite  material  in  order  to  obtain  an  energy  release  rate  per  unit  of 
newly  opened  crack  area.  It  is  this  figure  which  is  to  be  compared  with  critical  values 
for  a  particular  material. 

4.2.3  Programming 

The  foregoing  calculations  were  programmed  for  use  with  UAI/NASTRAN  for  analy¬ 
sis  and  I-DEAS  for  generation  and  display.  A  Fortran  program  was  written  to  perform 
the  following  functions: 

1.  Accept  user  definitions  of  the  Loner  and  outer  node  rings  as  “groups”  generated 
in  I-DEAS. 

2.  Obtain  displacements  for  these  nodes  from  a  “punch”  file  written  by  the  main 
NASTRAN  stress  analysis. 

3.  Open  a  new  NASTRAN  bulk  data  file  to  be  used  by  the  DMAP  program  de¬ 
scribed  below  in  calculating  energy  release  rates. 

4.  Write  displacements  to  this  file  in  DMIG  format. 

5.  Write  GRID  entries  for  all  the  ring  nodes,  CQUAD4  entries  for  the  wring  ele¬ 
ments,  and  associated  property  and  material  entries. 

6.  Write  executive  control  and  case  control  data  to  generate  thermal  loads  and  to 
access  the  required  DMAP  program. 

7.  Execute  the  NASTRAN  nm. 

8.  Find  and  display  the  strain  energy  release  rate  in  the  print  file. 

The  DMAP  program  consisted  of  about  30  lines  that  perform  the  following  oper¬ 
ations: 

1.  Read  node  and  element  data  and  generate  the  sti&ess  matrix  K  for  the  ring 
elements  only. 

2.  Generate  thermal  loads. 

3.  Read  displacements  d  from  DMIG  bulk  data  entries. 

4.  Compute  elastic  forces  by  the  matrix  miiltiplication  f  =  Kd. 

5.  Compute  a  thermal  load  vector  g  from  temperature  data,  thermal  expansion 
coefiScients,  and  element  geometry. 


4.3  Results 


Figure  24:  Intralayer  Embedded  Strain  Energy  Release  Rates 
6.  Compute  and  print  G  =  (f  +  g)^r. 

There  was  not  sufficient  time  in  Phase  I  to  demonstrate  virtual  crack  extension 
on  test  cases.  This  will  be  pursued  in  Phase  II. 

4.3  Results 

4.3.1  Intralayer  Embedded  Element 

Strain-energy  release  rates  predicted  by  the  intraJayer-embedment  model  are  shown 
in  Figure  24.  The  mode  I  energy  release  rate  is  considerably  higher  at  crack  initiation 
than  the  mode  II  release  rate.  However,  as  the  crack  grows  beyond  a  small  length, 
mode  II  release  rates  dominate. 

Note  that  the  energy  release  rates  shown  in  Figure  24  are  significantly  lower  than 
critical  strain  energy  release  rates  for  conventional  polymer  composites.  Typical  com¬ 
posite  interlaminar  strain  energ;;  release  rates  fall  in  the  range  250-450  N/m  (~1. 5-3.0 
Ib/in).  The  peak  rates  calculated  by  the  present  analyses  fall  some  20,000  times  be¬ 
low  that  range.  However,  the  intralayer-embedded  finite-element  model  considered 
applied  tensile  stresses  of  only  254  psi,  a  very  low  level  relative  to  the  composite’s 
capability.  Wisnom’s  analysis  shows  that  one  should  expect  mode  II  energy  release 
rates  to  be  proportional  to  the  square  of  the  applied  tensile  stress.  Thus,  critical 
energy  release  rates  mi^t  be  achieved  in  the  modeled  laminate  for  tensile  loadings  of 
the  order  of  38-45  ksi,  corresponding  to  roughly  half  the  strength  of  the  unidirectional 
composite  material. 


Mode  I  energy  release  rates  are  determined  chiefly  by  interlaminar  tension  rather 
than  shear.  Review  of  Table  2  reveals  that  stress  contributions  of  piezoelectric  and 
applied  tensile  loadings  are  approximately  equal  at  a  level  of  5000  psi  applied  axial 
stress.  Thus,  most  of  the  mode  I  energy  release  at  crack  initiation  is  due  to  the 
action  of  the  piezoelectric  element  for  the  relatively  low  applied  stress  of  254  psi. 
However,  the  mode  I  energy  release  rate  will  be  made  critical  only  by  increasing  the 
applied  loading,  since  the  induced  strain  loading  applied  in  the  model  is  very  close  to 
the  maximum  available  from  the  piezoelectric  material.  Assuming  that  mode  I  energy 
release  rate  scales  with  the  square  of  interlaminar  tensile  stress,  an  elementary  analysis 
shows  that  critical  energy  release  rates  might  be  achieved  for  applied  loadings  between 
11-15  ksi  (76-104  MPa).  Such  loadings  would  result  in  initiation  of  a  delaminating 
crack,  but  would  contribute  little  to  its  growth  beyond  a  short  length  because  of  the 
rapid  decay  of  mode  I  energy  release  rate,  as  shown  in  Figure  24. 

As  mentioned  above,  an  asymptotic  ^^ue  for  mode  II  energy  release  rate  can 
be  computed  based  on  the  applied  stress.  Figme  25  shows  the  ratio  of  the  mode 
II  energy  release  rate  to  its  asymptotic  value  for  the  intralayer-embedded  element 
case.  The  solid  line  in  this  figure  represents  the  data  also  presented  in  Figure  24. 
The  dashed  line  is  a  semi-log  extrapolation  of  the  trend,  indicating  an  estimate  of 
the  crack  length  at  which  the  asymptotic  value  will  be  achieved.  Once  the  crack 
grows  to  the  critical  length,  applied  stresses  sufficient  to  exceed  would  cause 
the  delamination  to  propagate  catastrophically.  This  length  is  projected  to  be,  at 
minimum,  roughly  0.01  in  for  the  intralayer-embedded  configuration  analyzed.  As 
this  is  a  small  value  relative  to  the  sensitivity  of  inspection  methods,  applied  loadings 
should  be  kept  below 

cTi  <  AGii^Eu  -  ■  (26) 

for  reliability.  Here,  h  is  the  total  thickness  of  the  laminate,  and  t  is  the  aggregate 
thickness  of  all  intelligent-material  elements  embedded  through  the  thickness  at  the 
critical  location.  For  the  case  analyzed  presently,  the  stress  limit  falls  in  the  range 
120-170  ksi,  or  roughly  a  factor  of  three  above  the  stress  level  at  which  delamination 
is  predicted  to  begin  to  be  governed  by  mode  II  behavior. 

4.3.2  Interlayer  Embedded  Element 

In  the  case  of  interlayer  embedment,  two  locations  in  the  intelligent,  composite  mar 
teriaJ  system  are  potentially  critical.  These  locations  are  shown  in  Figure  26,  and 
lie  along  the  interface  between  intelligent  and  host  composite  materials,  and  at  the 
point  where  the  host  composite  material  is  joined  beyond  the  extent  of  the  embedded 
element.  Strain-energy  release  rates  at  both  of  these  locations  were  investigated.  The 
model  considered  applied  axial  loadings  of  ±25  ksi. 

Mode  I  and  mode  II  strain-energy  release  rates  predicted  for  delamination  of  the 


Figure  25:  Intralayer-Embedded  Mode  II  Strain  Energy  Release  Rate  Ratio 
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Figure  26:  Interlayer-Embedded  Assumed  Crack  Locations 
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Delamination  Crack  Length  (in) 

Figure  27:  Strain  Energy  Release  Rates  along  Piezoelectric  Wafer  Interface 

embedded  element  from  the  host  composite  (loaded  in  compression)  are  shown  in 
Figure  27.  This  figure  shows  that  much  higher  values  are  obtained  in  mode  II  than 
mode  I,  although  the  model  predicted  these  values  somewhat  erratically.  Further 
investigation  revealed  that  the  oscillatory  predictions  were  due  to  oscillations  in  the 
calculated  values  of  relative  crack-sliding  displacements.  A  least-squares  fit  of  the 
predictions  is  also  illustrated  in  this  figure.  This  shows  an  onset  of  unstable,  mode  II 
crack  growth  up  to  a  length  of  roughly  0.001  in,  followed  by  stable  growth. 

Figure  28  shows  strain-energy  release  rates  calculated  at  the  point  where  composite 
layers  are  brought  together  outside  the  embedded-material  region.  In  the  present 
model,  this  location  lies  on  a  plane  of  enforced  symmetry,  so  crack  growth  is  possible 
only  mode  I. 

An  interesting  result  was  foimd  for  the  configuration  analyzed.  Crack  initiation 
must  occur  under  axial  compressive  loadings,  as  this  loading  produces  interlaminar 
tension  at  the  crack  tip.  However,  after  growing  a  small  distance,  the  stress  field  at 
the  crack  tip  reverses,  and  the  crack  grows  only  \mder  tensile  applied  loads.  Altho\igh 
figure  28  shows  this  subsequent  crack  growth  to  be  unstable,  the  linear  model  did  not 
consider  the  effects  of  crack-dosing  displacements.  As  axial  tension  is  applied  to 
the  configuration,  the  crack  closes  at  the  original  point  of  initiation,  lowering  tensile 
stress^  at  the  crack  tip.  This  ^ect  must  be  investigated  using  a  nonlinear  model  to 
characterize  strain-energy  release  rates  under  applied  tensile  loadings. 

Examination  of  figures  27  and  28  do  show  that  strain-energy  release  rates  for 
delamination  of  the  embedded,  intelligent  material  element  from  the  host  composite 
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Figure  28:  Strain  Energy  Release  Rates  along  Composite  Interface 


material  axe  much  higher  than  for  delamination  within  the  host  material  itself.  This 
delamination  is  dominated  by  mode  II  behavior,  and  would  be  expected  to  reach  Guc 
at  an  applied  stress  level  of  roughly  100  ksi. 


5.  Teusk  4  - 

Candidate  Experimental  Configurations 


5.1  Requirements 

The  integrity  of  intelligent  material  systems  in  service  will  \iltimately  be  ensured  by  a 
base  of  analytical  and  practical  experience  validated  by  extensive  testing.  The  objeo 
tives  of  any  testing  program  must  include  not  only  validation  of  theories  and  models, 
but  also  discovery  of  unknown  factors  influencing  material  performance.  While  this 
program  has  investigated  some  of  the  salient  aspects  of  intelligent  materials  integrity, 
a  subsequent  program  of  experimentation  should  be  undertaken  to  assess  the  validity 
of  the  present  results  (and  similar  approaches  pursued  by  others),  and  to  determine 
the  engineering  soxmdness  of  the  assumptions  made  in  such  models.  Accordingly,  this 
section  examines  some  basic  requirements  and  approaches  for  experimental  investi¬ 
gations  of  the  integrity  of  intelligent  materials  and  structures. 

First  and  foremost,  experimental  articles  should  be  built  to  replicate  configurations 
likely  to  be  employed  in  practice.  These  may  include  intralayer,  interlayer  or  fiber- 
embedded  elements,  surface-attached  elements,  or  even  elements  incorporated  directly 
into  the  host  composite  material  (e.g.,  whiskers).  Conventional  structural  composite 
materials  should  be  employed  as  host  materials. 

Second,  experimental  articles  should  be  subjected  to  the  environments  and  condi¬ 
tions  they  are  likely  to  see  in  service.  These  include  thermal,  atmospheric  and  loading 
environments  (both  static  and  dynamic). 

Third,  sensing  schemes  must  be  developed  to  track  the  occurrence  and  progression 
of  damage  within  the  intelligent  material  system  as  it  is  under  test.  Tracking  growth 
of  internal  damage  will  enable  detailed  validation  of  analytical  models,  and  will  bring 
to  light  any  important  factors  that  such  models  may  not  treat  well. 

Finally,  and  perhaps  in  the  longest  term,  the  experimental  program  should  build 
a  database  of  failure  modes  and  strengths  to  establish  reliable  engineering  properties 
for  intelligent  materials  and  structures. 

5.2  Sensing  Approaches 

A  key  issue  in  any  characterization  program  is  the  fidelity  of  sensing  elements  em¬ 
ployed  in  the  experiment.  The  idea  of  photoelastic  or  other  interferometric  sensing  of 
surface  strains  or  s\irface  “bulging”  suggests  itself.  However,  the  studies  performed 
so  fax  predict  very  small  surface  bulging  or  changes  in  smface  strains,  so  that  this 
approach  might  not  detect  a  crack  imtil  it  had  gotten  rather  far  along. 

Another  significant  finding,  as  discussed  above,  was  that  applied  loading  has  a 
greater  influence  than  piezoelectric  actuation.  This  suggest  that  procedures  “conven- 
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tionally”  used  for  detection  of  delamination  would  also  be  applicable  to  delaminations 
associated  with  embedded  elements. 

Through-the-thickness  acoustical  measurements  should  also  be  considered.  The 
rather  small  layer  thicknesses  suggest  possible  problems  with  spatial  resolution  in  this 
approach  (i.e.,  acoustic  wavelengths). 

In  Phase  II,  the  studies  initiated  in  Phase  I  will  be  continued,  partially  for  the 
purpose  of  designing  a  suitable  experiment.  Three-dimensional  models  may  be  used 
in  order  to  predict  two-dimensional  surface  responses.  Also,  we  will  research  current 
sensing  methodologies,  both  interferometric  and  acoustical. 


5.3  Discussion 

Perhaps  the  most  important  physical  quantity  to  be  assessed  by  experiments  in  this 
area  is  the  value  of  the  critical  strain-energy  release  rate  between  the  dissimilar  intel¬ 
ligent  and  host  composite  materials.  This  quantity  will  allow  assessment  of  loadings 
and  crack  dimensions  that  may  cause  catastrophic  propagation  of  interface  cracks 
which  decouple  the  intelligent  material  element  from  its  host  material.  Also,  the  pos¬ 
sibility  exists  that  the  critical  mode  II  energy  release  rate  for  a  candidate  ceramic 
active  material  may  dominate  the  interfacial  failure.  In  this  case,  the  ceramic  may 
be  found  to  have  lower  capability  in  this  mode  than  the  host  composite  material, 
causing  loads  to  be  limited  to  well  below  the  capacity  of  the  host  composite  in  order 
to  preserve  the  active  integrity  of  the  intelligent  material  system. 

If  energy-release  values  Gu^  are  foimd  to  be  lower  within  the  host  material  than 
across  host-intelligent  material  interfaces,  then  experimental  characterization  of  flaw 
propagation  due  to  external  loadings  need  only  consider  the  presence  of  a  stiff  element 
mimicking  the  intelligent  material  since  its  activation  plays  a  relatively  minor  role  in 
interlaminar  crack  propagation  in  the  host  material. 
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6.  Conclusions 


6.1  Observations 

The  following  observations  may  be  drawn  from  the  results  of  this  study: 

1.  The  presence  of  stiff,  active  material  elements  in  laminated  comp<^tes  causes 
significant  concentrations  of  interlaminar  tension  and  shear  stre^es  and  strains 
at  the  interfaces  between  the  two  materials.  The  magnitudes  of  tension  and 
shear  concentrations  caused  by  inducing  strains  in  the  active  element  are  roughly 
equivalent.  However,  the  magnitude  of  the  interlaminar  shear  stress  concentra¬ 
tion  is  roughly  twice  that  of  the  interlaminar  tension  stress  concentration  when 
in-plane,  external  loadings  are  applied  to  the  host  material. 

2.  For  cases  in  which  significant  external  loadings  are  applied,  local  interfacial 
stress  gradients  and  concentrations  are  dominated  by  the  magnitude  of  the 
external  loading.  The  amoxmt  of  strain  that  can  be  induced  in  conventional 
ceramic  actuator  materials  is  low  relative  to  the  capability  of  the  conventional 
composite  materials  (<  0.1%  vs.  ~  1%).  Hence,  strength-dominated  composites 
will  be  influenced  more  by  the  simple  presence  of  dissimilar,  intelligent  materials 
than  they  will  be  activation  of  these  materials.  Note  that  this  situation  may 
differ  for  intelligent  materials  capable  of  higher  activation  strains  (e.g.,  shape- 
memory  alloys). 

3.  Interfacial  stress  gradients  between  intelligent  and  host  composite  materials 
occur  over  distances  as  small  as  4-6  wafer-layer  thicknesses  from  the  interfaces. 

4.  For  intelligent  composites  dominated  by  external  loadings  (e.g.,  aircraft,  mis¬ 
siles,  launch  vehicles),  interfacial  stress  concentrations  are  minimized  ly  aligning 
the  composite  material  according  to  conventional,  strength-design  procedures. 
For  composites  dominated  by  intelligent-material  actuation  (e.g.,  satellite  struc¬ 
tures),  interfacial  shear  stress  concentrations  are  lowered  by  increasing  trans¬ 
verse  shear  stiffness  in  the  composite. 

5.  For  intelligent  material  wafers  embedded  within  composite  lamina  (intralayer 
embedment),  interface-crack  initiation  is  dominated  ly  mode  I  behavior,  but 
mode  I  stress  intensity  falls  off  rapidly  as  the  crack  advances.  Propagation  of 
interlaminar  failure  is  governed  by  model  behavior,  with  relatively  the  same 
characteristics  as  propagating  delaminations  in  conventional  composites.  Com¬ 
paring  mode  II  strain  energy  release  rates  to  their  asymptotic  values  allows  an 
estimate  of  the  critical  crack  length  to  be  made. 

6.  For  intelligent  matorial  wafers  embedded  between  layers  of  the  host,  compos¬ 
ite  material  (interlayer  embedment),  both  initiation  and  propagation  of  cracks 
along  the  material  interface  are  governed  by  mode  U  behavior.  Interlaminar 
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interfaces  near  the  embedded,  active  material  but  occurring  between  hc^t  ma- 
terial  layers  can  be  influenced  by  either  mode  I  or  mode  II  behavior,  although 
the  cases  analyzed  here  admitted  only  mode  I  crack  growth. 

7.  Both  tensile  and  compressive  loadings  are  important  for  crack  growth  in  inter¬ 
layer  embedded  material  systems.  Interlaminar  cracking  in  some  areas  may  be 
initiated  in  compression,  but  tend  only  to  propagate  in  tension. 

8.  Crack-closing  displacements  are  important  in  interlayer  embedded  conflgura- 
tions  due  to  the  coupling  of  in-plane  and  out-of-plane  forces  by  the  geometric 
configuration. 

9.  Strain-energy  release  rates  are  significantly  higher  along  the  interface  between 
intelligent  and  host-material  elements  that  between  host-material  layers  away 
from  the  embedded  element.  As  such,  critical  strain-energy  release  rates  along 
the  material  interface  must  be  characterized  if  the  integrity  of  this  interface  is 
to  be  assessed. 

6.2  Evaluation  of  Approach  and  Results 

The  approach  taken  in  the  present  study  has  illuminated  some  key  factors  influencing 
the  integrity  of  intelligent  materials  systems  and  structures.  In  particular,  that  in- 
terfacial  crack  propagation  is  governed  chiefly  by  interface  shearing  (mode  II)  action, 
and  that  interfacial  stress  concentrations  are  limited  to  material  regions  within  a  few 
thicknesses  of  an  embedded,  active-material  element. 

The  approach  taken  to  arrive  at  these  results  combined  analytical  and  finite- 
element  techniques  to  assess  responses  of  candidate  intelligent  material  configurar 
tions.  The  analytical  model  proved  appropriate  for  assessing  local  stress  gradients, 
and  admitted  rapid  analysis  of  various  material  configurations.  The  finite-element 
models  were  sufficiently  flexible  to  allow  interface  constraints  to  be  easily  modified 
for  assessment  of  strain-energy  release  rates.  However,  both  of  these  models  can  be 
improved  to  provide  more  useful  results. 

Two  specific  improvements  to  the  analytical  model  would  make  it  much  more 
useful  for  this  problem.  First,  the  assignment  of  Lagrangian  constraints  should  be 
made  at  the  modeling  level,  rather  than  at  the  numerical  level  as  is  now  the  case. 
This  would  allow  engineers  to  specify  particular  connectivity  topologies  of  interest, 
rather  than  just  the  fully  connected,  rectilinear  topology  implemented  to  date.  The 
analytical  formulation  would  then  become  useful  for  studies  of  stress  distributions 
around  interfacial  cracks.  Second,  higher-order,  nonlinear  strain-displacement  should 
be  implemented  to  assess  the  effects  of  geometric  nonlinearities  on  local  stresses  near 
interfaces.  The  polynomial  natme  of  the  assumed  displacement  (and  strain)  functions 
in  the  analytical  solution  lend  themselves  easily  to  sudi  atigmentation. 


The  linear  finite-element  models  used  in  the  present  study  should  be  modified 
to  include  nonlinear  effects  associated  with  geometric  stiffness  (particularly  for  the 
interlayer-embedded  case)  and  with  gaps  resulting  from  crack  opening  displacements. 
Modeling  extensive  flaw  propagation  in  the  interlayer  embedded  configurations  under 
cyclic  loading  requires  gap  constraints  to  be  considered  since  extension  of  cracks  near 
the  transition  region  are  associated  with  opening  and  closing  displacements.  Accurate 
calculation  of  strain-energy  release  rates  in  these  cases  requires  a  kinematically  correct 
model  of  deformations  between  internal,  free  surfaces. 

Finally,  if  the  present  techniques  are  to  provide  guidelines  for  engineering  design  of 
intelligent  structures,  experimental  data  is  needed  to  assess  mode  II  strain-energy  re¬ 
lease  rates  along  interfaces  between  conventional,  structural  composite  materials  and 
practical  intelligent  materials  (e.g.,  piezoceramics,  shape-memory  alloys).  Without 
these  data,  analytical  predictions  will  lack  the  accuracy  needed  for  their  successful 
application. 


